Depairing critical current achieved in superconducting thin films with 
through-thickness arrays of artificial pinning centers 
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Large area arrays of through-thickness nanoscale pores have been milled into superconducting Nb 
thin films via a process utilizing anodized aluminum oxide thin film templates. These pores act as 
artificial flux pinning centers, increasing the superconducting critical current, J c , of the Nb films. By 
optimizing the process conditions including anodization time, pore size and milling time, J c values 
approaching and in some cases matching the Ginzburg-Landau depairing current of 30 MA/cm 2 at 
5 K have been achieved - a J c enhancement over as-deposited films of more than 50 times. In the 
field dependence of J c , a matching field corresponding to the areal pore density has also been clearly 
observed. The effect of back-filling the pores with magnetic material has then been investigated. 
While back-filling with Co has been successfully achieved, the effect of the magnetic material on J c 
has been found to be largely detrimental compared to voids, although a distinct influence of the 
magnetic material in producing a hysteretic J c versus applied field behavior has been observed. This 
behavior has been tested for compatibility with currently proposed models of magnetic pinning and 
found to be most closely explained by a model describing the magnetic attraction between the flux 
vortices and the magnetic inclusions. 



J 300 nm V* , •% \**A 




I. INTRODUCTION 

In recent years, the state-of-the-art critical current 
density, J c , of superconducting films has improved sig- 
nificantly as researchers have developed nanocompos- 
ites incorporating artificial pinning centers pQ. Non- 
superconducting inclusions, such as insulating particles 
within a YBa 2 Cu 3 07_5 matrix [2, 3 , serve to pin mag- 
netic flux. Typically, such inclusions are incorporated 
randomly into the superconducting matrix, without the 
fine control of size and spacing that is necessary for op- 
timization of J c . The first part of this work describes 
a versatile method for fabricating artificial pinning cen- 
ters with a highly controlled geometry: nanoscale pores, 
shown in Fig. [T] are transferred from an anodized alu- 
minum oxide (AAO) template into a superconducting 
film. Anodization of alumina is a well-studied and re- 
liable self-assembly process [2HZ] that can create billions 
of nanoscale pores in parallel over a large area, with con- 
trol over both the spacing and size of the pores. Such a 
scalable process could be integrated into the commercial 
production of long lengths of superconducting wire [I], 
unlike slow, serial patterning processes such as electron 
beam lithography and focused ion beam patterning. 

The size and spacing of AAO pores is well matched 
to the vortex lattice that forms in practical type-II su- 
perconductors. The spacing can be varied from at least 
300 nm to below 50 nm, corresponding to the spacing of 
the vortex lattice at fields from 30 mT to 1 T. There have 
been several efforts to exploit this natural match between 
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FIG. 1. (Color online) Anodized aluminum oxide (AAO) 
forms nanoscale pores over large areas, as shown in this scan- 
ning electron micrograph of the surface of an AAO film on 
Nb (schematic, upper inset). The image shows a disordered, 
natural arrangement of pores of well defined diameter (75 ± 
10 nm). The rings in the two-dimensional Fourier transform 
of the image (lower inset) indicate a lack of crystalline order, 
but a defined pore spacing (~140 nm). 



AAO and the vortex lattice, including depositing super- 
conducting niobium (Nb) on top of AAO 8 , depositing 
NbN on wires templated by AAO [9], and using AAO as 
a mask to ion mill holes in Nb [TD]. Here, we further de- 
velop this last process, creating very strong pinning cen- 
ters by ensuring through-thickness holes in the Nb layer, 
and improving the process reliability through secondary 
ion mass spectrometry during milling. 
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Having created strong and well-characterized pinning 
centers in the form of holes in Nb, we then address 
the basic scientific question of whether magnetic pin- 
ning centres can improve on their non-magnetic equiv- 
alents, a source of significant debate given the impor- 
tance of progressing beyond the capabilities presently of- 
fered by core pinning approaches The AAO sys- 
tem is ideal for backfilling the holes with magnetic ma- 
terial, leading to through-thickness magnetic inclusions 
that can be directly compared to their non-magnetic 
counterpart. While various mechanisms have been pro- 
posed for magnetic pinning [T^HTB] . the most technolog- 
ically promising |17j requires such through-thickness in- 
clusions, rather than magnets placed above or below a 
superconducting film, as in most of the experiments con- 
ducted thus far [TM2"T] . In previous experiments with 
through-thickness magnetic inclusions [22H24] , such in- 
clusions inevitably affect (and usually increase) the core 
pinning of vortices in the film, making it impossible to 
isolate the contribution of magnetic interactions to J c . 
In the system presented here, we overcome this difficulty 
by changing the magnetization of through-thickness in- 
clusions in situ and observing the resultant change in J c . 
This relationship between pinning and magnetization is 
manifested as a novel, hysteretic behavior of J c under an 
applied field. 



II. NANOFABRICATION OF ARTIFICIAL 
PINNING SITES 

The cross-sectional structure of a typical film is illus- 
trated schematically in Fig. [2]ja). A 36 nm thick su- 
perconducting niobium (Nb) layer and a 670 nm thick 
aluminum (Al) layer are deposited in a single sputtering 
run by moving the oxidized silicon substrate past two DC 
magnetrons. The chamber is baked out before deposition 
and liquid nitrogen cooled during deposition, resulting in 
a background pressure of less than 10 -8 mbar. The alu- 
minum layer is then anodized into AI2O3 in 0.3 M oxalic 
acid at room temperature following a procedure for form- 
ing porous AAO on silicon [6]. During anodization, the 
substrate and film edges are protected from the acid bath 
with a temporary adhesive (glycol phthalate). 

As shown in Fig. [5Jb), the anodization process is self- 
terminating. A constant voltage is applied, and the cur- 
rent indicates the rate of anodization. The anodization 
rate is initially high but drops rapidly as the entire top 
surface anodizes to form AI2O3. There is a resurgence in 
the anodization rate as pores begin to form and the ex- 
posed surface area increases. Once all pores are formed, 
the process reaches a roughly steady state and the an- 
odization front progresses through the thickness of the 
film. The anodization rate drops sharply once the an- 
odization front reaches the bottom of the Al layer, al- 
though anodization is continued until the rate begins to 
level off again, to allow the pores to reach close to the bot- 
tom of the Al layer, separated only by a residual barrier 



layer which is not removed even under further anodiza- 
tion. By varying the anodization voltage, variable pore 
spacings can be achieved. In this work, voltages of 40 
and 57 V were used, resulting in average pore spacings of 
approximately 100 and 140 nm. After anodization, the 
pores are just 20 nm in diameter, too small for subse- 
quent pattern transfer by Ar + ion milling [7 . To widen 
the pores, we isotropically wet etch using phosphoric acid 
(5% aqueous solution) for 56, 75 and 83 minutes at room 
temperature, obtaining pores of diameter 60, 75 and 80 
±10 nm. A scanning electron microscope (SEM) image 
of the cleaved edge of the structure (Fig.|3ja)) shows that 
after anodization and pore widening, the pores terminate 
about 60 nm above the Nb layer. 

To extend the pores through the Nb, we require an 
anisotropic etch (an isotropic wet etch would completely 
remove the Nb). We initially tried a CF4 plasma gener- 
ated by an RF field applied perpendicular to the sample 
surface (using a parallel plate capacitor geometry to ap- 
ply the field). This etch would be anisotropic and also 
chemically selective, etching Nb while not significantly 
etching AI2O3. SEM images indicated, though, that the 
Nb was not etched by this process. This may be due to 
the residual AI2O3 at the bottom of the pores (in spite 
of the phosphoric acid etch) preventing the plasma from 
reaching the Nb. 

We are, however, able to extend the pores using argon 
ion milling, the effectiveness of which has been demon- 
strated by milling pores into silicon substrates through 
AAO [7]. The cross-section of our film after milling is 
shown in Fig.[3jb). This technique etches the AAO, and 
in fact etches the surface more quickly than the bottom 
of the pores, presumably due to increased likelihood of 
redeposition of material inside the pores. Thus it is nec- 
essary to etch the template until the pores become suffi- 
ciently shallow that material is efficiently removed from 
them. SEM images indicate that the pores do not ex- 
tend fully through the Nb until the template is thinned 
to less than 200 nm. Using a beam voltage of 1000 V with 
an ion current density of 3.9 mA/cm 2 (20 mA nominal 
beam current), we mill for approximately 16 minutes to 
remove 820 nm of template, leaving 170 nm of template 
intact. The pores clearly extend through the Nb layer, 
and when the Nb is peeled away from the Si02 layer of 
the substrate, as in Fig. [3](c), the pores are seen to pen- 
etrate through into this substrate layer. We speculate 
that ion milling followed by CF4 etching could be used 
to produce through-thickness pores in much thicker films. 

Milling for a fixed time with nominally identical beam 
parameters was not sufficiently reproducible to reliably 
remove the desired amount of template without damag- 
ing the entire Nb layer. We overcame this problem using 
secondary ion mass spectrometry to monitor the sputter 
products during ion milling, as shown in Fig. [4] We were 
able to detect a downward kink in the milling rate of the 
Al from the template, accompanied by a slight increase 
in the Nb milling rate. By examining the cross section 
of milled samples with SEM, we established an empir- 
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FIG. 2. (Color online) (a) Schematic of the anodization process. The film's cross section (not to scale) is depicted immersed 
in an oxalic acid bath at a fixed potential, forcing oxidation of the Al and pore formation. The current measured during this 
process, shown in (b), reveals the various phases of pore formation, which self-terminates at the Nb layer. 



ical criterion of a 25% reduction in Al milling rate for 
terminating the milling with approximately 170 nm of 
template remaining. The increase in Nb rate during the 
last 1.7 minutes (90 nm of template milling) corroborates 
the idea that the template must be thinned to less than 
200 nm to effectively mill the Nb. In retrospect, then, 
the 990 nm template is unnecessarily thick. 



III. HIGHLY IMPROVED J c WITH TUNABLE 
MATCHING FIELD 

Milling the Nb layer through the AAO template results 
in a dramatic increase in critical current. The results 
shown in Fig. [5] are obtained from transport measure- 
ments on a 12 /im wide by 10 /im long bridge patterned 
by focused ion beam. The magnetic field H is applied 
perpendicular to the film surface. At 5 K, the milled 
sample's J c is more than 50 times higher than that of a 
non-anodized control sample. This calculation of J c does 
not even take into account the 30% reduction in film cross 
section resulting from the introduction of the pores. Thus 
the artificial pins created by the pores dominate over the 
natural pinning sites present in an as-deposited film. 

This J c enhancement cannot be dismissed as merely 
poor performance of the control sample. Multiple control 
samples from the same deposition run were compared and 
their J c values found to differ by at most 30%, typically 
5%. Between deposition runs, J c values differed by up to 
a factor of three, but these differences can be accounted 
for by corresponding variations in T c , which is known 
to be extremely sensitive to the background pressure of 
oxygen during deposition. 

What is more remarkable is that at 5 K (T c =7.5 K), 



the critical current at zero field reaches the calculated 
Ginzburg-Landau depairing current [25], 
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where fio = 4ir x 10~ 7 H/m is the permeability of free 
space, $ = h/2e = 2.07 x 10~ 15 Wb is the flux quan- 
tum, and the upper critical field H C 2 is taken from the 
data as the field at which the IV response becomes linear 
(Ohmic), corresponding to the transition into the nor- 
mal state. The penetration depth at K for a film of 
this thickness in accordance with dirty limit BCS theory, 
A(0) = 120 nm is taken from Ref. and scaled for 
temperature according to the two-fluid model [25], 
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This astounding pinning performance is only expected 
when each vortex sits on an ideal, columnar hole. Al- 
though we did not obtain this performance consistently 
across samples (the best sample is shown), it is impres- 
sive that the sample approaches this limit at all. Such 
high performance only appears in one recent work, also 
on nanostructured films [2"7] . 

J c shows significant additional enhancement around 
0.10 ±0.01 T and multiples thereof (most evident at 7 K 
or on the linear scale of Fig. [8]). This matching field rep- 
resents a filling of one vortex per pore, which accords 
with the pore density calculated from SEM images of the 
film surface, counting pores over a known area. 
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FIG. 3. (Color online) SEM images illustrate the pore transfer 
process, (a) An image of the cleaved edge of the film before 
milling, (b) After milling, the AI2O3 layer is thinned and 
the pores extend through the Nb layer, (c) Peeling away the 
Nb and AI2O3 layers reveals pits in the SiC>2 substrate layer, 
proving that the pores extend through the Nb. 



The matching field is unchanged between 5 and 6 K, but 
shifts significantly to lower fields at the highest measure- 
ment temperature, 7 K. At this temperature, the coher- 
ence length, £ (calculated as £ = \f<&o/H C 2) increases to 
90 nm (from 40 and 50 nm at 5 and 6 K, respectively), be- 
coming comparable to the pore spacing of about 145 nm. 
The more closely spaced pairs of pores may then begin to 
act as single pinning sites, decreasing the effective den- 
sity of pinning sites and thereby lowering the effective 
matching field. 

The matching field can be manipulated by changing 
the pore density, which is easily controlled via the linear 
relationship between pore spacing and anodization volt- 
age [B] . Fig. [6] compares results obtained from transport 
measurements on bridges 120 /im wide by 500 /jm long 
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FIG. 4. (Color online) Secondary ion mass spectrometry pro- 
vides a reliable trigger for terminating the milling process. 
Once the AI2O3 layer is sufficiently thin, the Al rate de- 
creases rapidly, accompanied by an increase in the Nb rate 
as expected. 
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FIG. 5. (Color online) Milling pores (75 nm diameter, 140 nm 
spacing) in a Nb film dramatically improves its transport crit- 
ical current, compared with an unanodized control. Measure- 
ments were taken at 5 K, 6 K and 7 K (T c = 7.5 K). Within 
each set of curves, the thin solid line is the calculated depair- 
ing critical current, Jd, the thick solid line with points is the 
porous film, and the dashed line with points is the control. 



photolithographically patterned into films anodized at 40 
and 57 V. In both cases, the pore density calculated from 
the SEM images is consistent with the observed matching 
field. 
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FIG. 6. (Color online) Films anodized at different voltages 
have different pore densities (insets) which in turn result in 
different matching fields in critical current versus magnetic 
field. (57 V: 75 nm diameter, 140 nm spacing; 40 V: 50 nm di- 
ameter, 100 nm spacing.) The measurement temperature was 
6 K; T c = 7.2 K, 8.0 K, and 8.0 K for the 57 V, 40 V, and 
control films, respectively. J c values of the 57 V sample have 
been plotted reduced by a factor of two for ease of comparison 
with the other data. 



is therefore a result of the smaller-than-expected depo- 
sition area, and could be corrected by terminating the 
deposition sooner. However, we found no way of over- 
coming the non-uniform growth. 

We therefore chose, ultimately, to backfill the pores 
by electron beam evaporation of Co. Evaporation of a 
130 nm thick Co layer was performed in high vacuum, 
leading to directional deposition with low sidewall cov- 
erage. This is desirable in order to avoid clogging the 
pores before they are completely filled. Cross sectional 
SEM images, as in Fig. ^b), confirm that evaporation 
is effective, showing Co (which displays a grain structure 
that distinguishes it from the A AO and Si02) present 
in the bottoms of the pores, penetrating through the 
Nb layer. As expected, Co is also deposited on top of 
the template. However, we are able to remove this top 
layer without affecting the inclusions by off-axis Ar + ion 
milling. Milling at 70° from the normal to the substrate 
with a beam voltage of 1000 V and an ion current density 
of 3.9 mA/cm 2 for five minutes removes the top Co layer 
and half of the remaining 170 nm of template. A cross 
section of the final film is shown in Fig. ff\c). 



V. EFFECT OF MAGNETIZATION ON 
CRITICAL CURRENT 



IV. FABRICATION OF MAGNETIC PINS 

We now use this system of strong artificial pinning cen- 
ters to test whether magnetic inclusions can improve on 
the existing core pinning, as discussed in the Introduc- 
tion. Backfilling the pores with magnetic material pro- 
vides a comparison between magnetic and non-magnetic 
pins, without changing the superconducting material or 
geometry of the pins. We chose to backfill with cobalt, a 
magnetic material that is straightforward to deposit by 
several methods, and has a large saturation magnetiza- 
tion, producing a flux density of 1.8 T. 

Initially, we electrodeposited Co from a cobalt sulfate 
solution (CoS0 4 -5H 2 1.3 M, H 3 B0 3 0.7 M, pH 2.0) at 
room temperature [2S]. Electrodeposition appears to be 
the ideal method of backfilling the pores, in that the Nb 
acts as an electrode while the insulating AAO template 
remains uncoated, so that the film is deposited only at 
the pinning sites as desired, rather than as a continuous 
sheet. As illustrated in Fig. ^a), however, the deposi- 
tion proved to be very non-uniform, growing in only a 
third of the pores. Furthermore, many of the Co growths 
overgrew the template after only 5 seconds of deposition. 
The deposition was terminated based on a thickness cri- 
terion of 150 nm, which translates to a charge (charge is 
the product of volume, atomic density of Co, and charge 
per atom); the charge is monitored as the integral of 
the current during deposition. This termination crite- 
rion fixes the deposition volume, so the thickness will in- 
crease as the deposition area decreases. The overgrowth 



Contrary to our initial expectations, adding Co to the 
pores does not significantly increase J c . To obtain the 
data in Fig. |8j films were coated with Co over part of 
the substrate using a shadow mask during evaporation. 
Transport bridges 120 /im wide by 500 /im long were 
then photolithographically patterned in each section of 
the sample, allowing a comparison of magnetic and non- 
magnetic J c on the same sample. For a film containing 
75 nm diameter pores (Fig. [8^ a)), J c decreases across all 
applied fields with the addition of Co. T c is 7.2 K for 
both of the bridges, and therefore does not offer an ex- 
planation for the decrease. Increasing the pore diameter 
is expected to decrease the overall J c as the supercon- 
ducting cross-section is reduced. In itself, this is insuffi- 
cient to explain the factor of three reduction in J c values 
on increasing the pore size of the non-magnetic samples. 
It could, however, explain the slight reduction in self- 
field J c when comparing the two magnetic samples. The 
magnetic sample with 80 nm diameter pores shows more 
closely comparable values to the non-magnetic sample, 
with even some enhancement at |/io-ff| > 0.3 T, but still 
a lower value at zero field, where Lorentz force reduction 
predicts the greatest increase |17j . The optimal pin size 
for a vortex core pin is well known to be ~2£ = 76 nm for 
Nb - the size at which the non-magnetic sample exhibits 
the greatest enhancement. The optimal pin size for mag- 
netic pinning depends on the mechanism at work. For 
Lorentz force reduction, the optimal pin size would be 
one that could carry a full quantum of flux. For fully 
saturated Co (/i M s = 1.8 T), this would be a 38 nm 
diameter. Based on the actual Co magnetisation loop 
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FIG. 7. (Color online) To test the effect of magnetization on 
pinning, we backfill pores with Co. (a) Electrodeposition of 
Co is non-uniform, filling some pores and not others, as seen 
in this SEM image, (b) Evaporation results in uniform filling, 
but also a layer of Co on top of the AI2O3. (c) The top layer 
of Co is removed by milling at a shallow angle, leaving Co in 
the pores. 

of these samples, however, an applied field of just above 
0.1 T [hqM ps 0.5 T) would be optimal for these pore 
diameters. At this field, both non-magnetic samples are 
outperforming their magnetic counterparts, contradict- 
ing the predictions of the Lorentz force reduction theory 
of magnetic pinning. 

There are likely factors other than the magnetism of 
the Co behind the seemingly disappointing performance 
of the Co-filled samples compared to their unfilled coun- 
terparts. For example, the core pinning potential is sig- 
nificantly reduced by the proximity effect when a conduc- 
tor, whether magnetic or not, is brought into contact with 
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FIG. 8. (Color online) A film backfilled with Co shows lower 
J c than its unfilled counterpart, possibly due to an adverse 
proximity effect. Two films are shown with different pore di- 
ameters (75 and 80 nm) but the same pore spacing (140 nm), 
measured at 6 K. An unanodized control sample is also shown. 
The T c of the control sample is 8.2 K; all others are 7.2 K. 



the superconductor (i.e. the Co deposited in the pores). 
The gradient in potential energy (and hence core pinning 
force) would then be determined by the length scale of 
this proximity effect's suppression of superconductivity, 
rather than the coherence length. We suggest, then, that 
depositing a ferromagnetic insulator would allow a more 
perfect experimental comparison than our ferromagnetic 
metal. 

Despite the overall decrease in J c , there is evidence for 
a magnetic contribution to pinning in the form of a hys- 
teresis in J c versus applied field that is not present in the 
non-magnetic samples (Fig. [9][a)). This hysteresis arises 
from the hysteretic magnetization curve of the Co inclu- 
sions. Fig. [9^b) shows the magnetization (at 300 K) of 
an identical sample that is not patterned into transport 
bridges. We expect that this magnetization is not sig- 
nificantly different in the presence of superconductivity, 
since the large demagnetization factor of a thin film in 
perpendicular field H makes the field inside and outside 
the sample nearly the same. Furthermore, the internal 
field is approximately homogeneous when the vortices are 
spaced more closely than the thin film effective penetra- 
tion depth (true for \fioff\ > 2 mT). 

The hysteresis in J c tells us that pinning is stronger 
when the inclusion magnetization is aligned with the vor- 
tices: for positive applied fields, the critical current is 
higher when descending in applied field (the blue curve 
in Fig. |9^a)), which corresponds to more positive mag- 
netization (the blue curve in Fig. gb)). The opposite 
is true for negative applied fields. This is qualitatively 
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FIG. 9. (Color online) (a) Hysteresis in J c versus applied field 
at 5 K indicates that aligning the magnetization of the Co 
inclusions parallel to the superconducting vortices increases 
J c . The hysteresis follows the hysteretic magnetization of the 
Co, shown in (b). A sample of Co inclusions with the top 
layer of Co removed is measured at 300 K with the applied 
field perpendicular to the film surface. Red and blue curves 
indicate increasing and decreasing applied field, respectively. 



consistent with the Lorentz force reduction theory [T7] . 
in which the force acting on a vortex pinned by a mag- 
netic inclusion is diminished by simultaneously acting on 
the inclusion, but is also consistent with several other 
mechanisms of magnetic pinning, including field compen- 
sation [14j whereby the applied magnetic field is to some 
extent compensated by the field generated by the mag- 
netic material, and magnetic attraction |15j in which the 
pinning force acting on the vortex is supplemented by an 
additional magnetic component. 

To distinguish between these theoretical mechanisms, 
we must look quantitatively at the size of the hysteresis. 
In Fig. [lOj we have re-plotted the hysteretic J c versus 
H for H > 0. We first attempt to fit the Lorentz force 
reduction theory to these data using Eqs. 8 and 10 of 
Ref. [T2]: 




'o 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
applied field (T) 



FIG. 10. (Color online) The size of the hysteresis in J c (data 
shown as a solid line with points; red curves correspond to 
increasing applied field, and blue to decreasing applied field) 
allows comparison with the calculated predictions of several 
theories of magnetic pinning. The large difference between 
J c on ascending and descending sweeps of the applied field 
expected from Lorentz force attraction (dashed curves) is not 
observed, nor is the substantial field shift predicted by the 
field compensation theory (dotted curves). The data is most 
consistent with the simple magnetic attraction calculation 
(dash-dot curves). 

B = B' + tx M(B') ■ f (5) 

where J c ' / J c is the ratio between the magnetic and non- 
magnetic critical current densities, / is the volume frac- 
tion of the inclusions (equal to the area fraction of the 
pores, measured from SEM images to be 0.16) [29 , B is 
the applied magnetic field, and B' is an effective field. 
This effective field is calculated by numerically inverting 
the second equation, using the measured magnetization 
curve of the inclusions, M(H). Having calculated the ra- 
tio J c '/J c as a function of H , we multiply it by the aver- 
age of the measured ascending and descending J c curves 
as an approximation of the form of the non-magnetic J c . 
(As explained above, the non-cobalt filled sample proved 
unsuitable for direct comparison with this sample.) The 
prediction of the theory in this case is not the absolute J c , 
since we have no direct measurement of the non-magnetic 
J c , but it does yield the ratio between ascending and de- 
scending curves. The measured data cannot satisfy the 
large predicted divergence of the ascending and descend- 
ing curves, however. In short, Lorentz force reduction 
theory predicts a much larger effect of magnetization on 
J c than that observed. 

The equivalent prediction for the field compensation 
scenario |14j requires transforming the field axis, H — > 
H' , according to the measured magnetization of the in- 
clusions. 



— = (1-/)- 
J c 1 ' B 1 



(4) 



H' = H- fM(H) (6) 
where, as before, / is the area fraction of the inclusions. 



The second term represents the spatially averaged re- 
turn field of the inclusions, which is subtracted from the 
applied field. As above, we employ the average of the 
ascending and descending data as the starting point for 
the transformation. The result is plotted as the solid 
lines in Fig. [lOj The hysteresis in the calculated curves 
follows from the hysteresis in the measured M(H). The 
large shift in field predicted by this model is also not 
consistent with the observed J c behavior. 

Finally, we consider the magnetic attraction between 
vortices and magnetized inclusions [15) . The pinning po- 
tential per length of vortex is modified by the addition of 
the Zeeman energy term —AMH/2, where A is the area 
of the inclusion perpendicular to the applied field direc- 
tion, M is its magnetization, and H is the field applied to 
the inclusion by the vortex lattice. The additional mag- 
netic pinning force is then the gradient of this additional 
pinning potential. Both M and H vary as a function of 
H, the macroscopic applied field. We took M(H) to be 
as shown in Fig. [9^b) . To estimate H, we summed the 
field from two neighboring vortices at a point between 
them a distance r from one of them: 



H{r) = 



27T/i A 2 



K„ 



A 



(7) 



where Kq is a modified Bessel function, and a is the ex- 
pected vortex lattice spacing for the particular applied 
field, 



2$n 



Mo 



Hs/2, 



(8) 



We then calculated the magnetic pinning force, /™ as , 
as the gradient of the resulting potential, evaluating the 
result with the inclusion offset from the center of the 
vortex by a distance £, as this is the position we expect 
a pinned vortex to occupy due to the strong core pinning 
of the inclusions: 



Jp 



mag 



AM 



dH 

dr 



(9) 



Below the first matching field (i.e. at vortex densities 
lower than one vortex per pore), we assume each vortex 
is paired with one inclusion, thus the additional pinning 
force above contributes directly to a change in critical 
current, AJ C = fp la9 /$o- Above this field, we applied 
the additional magnetic pinning force to the vortex lat- 
tice as a whole, dividing the single vortex pinning force 
by the ratio of the density of vortices to the density of 
pores. 

To obtain the magnetic attraction (dash-dot) curves 
shown in Fig. |10[ we add the calculated magnetic pinning 
force (which is positive when M and H are aligned, and 
negative when they are anti-aligned) to the average of 
the ascending and descending sweeps from the measured 
data. Thus we are comparing the difference between as- 
cending and descending J c curves in the data and the 



calculation. They are in fact comparable, particularly 
above the first matching field. Below the first matching 
field, the calculation yields a larger magnetic contribu- 
tion to pinning than that observed. Idealizations in the 
calculation such as the placement of the vortices entirely 
on inclusions may account for the disparity. 



VI. CONCLUSIONS 

We have exploited the scalable, self-assembled 
nanofabrication technique of porous anodized alumina to 
make through-thickness pores in a superconducting film. 
These pores act as highly effective artificial pinning cen- 
ters with a tunable size and density. With this system, we 
are able to unambiguously demonstrate and quantify the 
effect of magnetization on pinning, observing a unique 
hysteretic J C {H) characteristic stemming from the hys- 
teresis of the inclusion magnetization. We find, however, 
that the size of this effect is smaller than predicted by 
the Lorentz force reduction theory. This is not unex- 
pected, primarily because the theory is computed for an 
ideally soft ferromagnet whose magnetization is propor- 
tional to applied field whereas the Co nanopillars have 
sufficient coercivity to prevent their magnetization from 
being realigned by individual vortices. A field compen- 
sation theory similarly fails to describe the observed be- 
havior. Instead, our magnetic pinning is best modeled by 
considering the magnetic attraction between the vortex 
lattice and the inclusion lattice. 
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